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Abstract 

We search for the strange pentaquark 0+ using kaon interactions in the material of the Belle 
detector. No signal is observed in the pKg final state, while in the pK^ final state we observe 
~ 1.6 • 10^ A(1520) — > pK~ decays. We set an upper limit on the ratio of O"*" to A(1520) yields 
(t(G"'")/(t(A(1520)) < 2% at 90% CL, assuming that the 0"*" is narrow. We also report on searches 
for strange and charmed pentaquarks in B meson decays. These results are obtained from a 
155 fb~^ data sample collected with the Belle detector near the T(4S') resonance, at the KEKB 
asymmetric energy e^e" collider. 

PACS numbers: 13.25.Hw, 13.75.Jz, 14.20.Jn, 14.20.-c 
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INTRODUCTION 



Until recently, all reported particles could be understood as bound states of three quarks 
or a quark and an antiquark. QCD predicts also more complicated configurations such as 
glueballs gg, molecules qqqq and pentaquarks qqqqq. Recently, observations of the pen- 
taquark B"*" = uudds have been reported in the decay channels K~^n P and pKs 0. 
Many experimental groups have confirmed this observation and the isospin 3/2 members of 
the same pentaquark multiplet have also been observed j^. Evidence for the charmed pen- 
taquark 0° = uuddc has also been seen j^. The topic attracts enormous theoretical interest. 
However the existence and properties of pentaquarks remain a mystery. Some experimen- 
tal groups do not see the pentaquark signals. The non-observing experiments correspond 
to higher center-of-mass energies. It has been argued that pentaquark production is 
suppressed in the fragmentation regime at high energies. 

Charged and neutral kaons are copiously produced at Belle. We treat kaons as projectiles 
and the detector material as a target, and search for strange pentaquark formation, KN — >■ 
G"*", and production, KN — >■ G^X. The kaon momentum spectrum is soft, with a most 
probable momentum of only 0.6GeV/c. Therefore we can search for 0^ formation in the 
low energy region. 

We also search for strange and charmed pentaquarks in the decays of B mesons, where the 
suppression of production observed in s channel e+e~ collisions p] may be absent. Studies 
of B meson decays have proved to be very useful for discoveries of new particles (such as 
P-wave cq states), therefore it is interesting to search for pentaquarks in B decays although 
no firm theoretical predictions for branching fractions exist. 

DETECTOR AND DATA SET 

These studies are performed using a data sample of 140 fb^^ collected at the T(45') 
resonance and 15fb~^ at an energy 60MeV below the resonance. The data were collected 
with the Belle detector ^ at the KEKB asymmetric energy e^e~ storage rings |^. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a three 
layer silicon vertex detector (SVD), a 50-layer cylindrical drift chamber (CDC), a mosaic of 
aerogel threshold Cherenkov counters (ACC), a barrel-like array of time-of-flight scintillator 
counters (TOF), and an array of CsI(Tl) crystals (ECL) located inside a superconducting 
solenoidal coil that produces a 1.5 T magnetic field. An iron flux return located outside the 
coil is instrumented to detect muons and mesons (KLM). 

The proton, kaon and charged pion candidates are identified based on the dE/dx, TOF 
and Cherenkov light yield information for each track. Ks candidates are reconstructed via 
the vr+TT^ decays and must have an invariant mass consistent with the nominal Ks mass. 
The Ks candidate is further required to have a displaced vertex and a momentum direction 
consistent with the direction from its production to decay vertices. 

SEARCH FOR 6+ USING KAON INTERACTIONS IN THE DETECTOR MATE- 
RIAL. 

The analysis is performed by selecting pK~ , pK^ and pKs secondary vertices. The 
protons and kaons are required not to originate from the region around the run-averaged 
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FIG. 1: The xy distribution of secondary pK vertices for the barrel (left) and endcap (right) 
parts of the detector. 

interaction point (IP). The proton and kaon candidate are combined and the pK vertex is 
fitted. The xy distribution of the secondary pK~ vertices is shown in Fig. ^ for the barrel 
part (left) and for the endcap part (right) of the detector. The double wall beam pipe, three 
layers of SVD, the SVD cover and the two support cylinders of the CDC are clearly visible. 
The xy distributions for secondary pK^ and pKg vertices are similar. 

The mass spectra for pK~ , pK^ and pKs secondary vertices are shown in Fig. |21 No 
significant structures are observed in the M{pK~^) or M{pKs) spectra, while in the M{pK~) 
spectrum a A(1520) signal is clearly visible. We fit the pK~ mass spectrum to a sum of a 




FIG. 2: Mass spectra of pK~^ (left), pK (right, points with error bars) and pKs (right, histogram) 
secondary pairs. The fit is described in the text. 

A(1520) probability density function (p.d.f.) and a threshold function. The signal p.d.f. is a 
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D-wave Breit-Wigner shape convolved with a detector resolution function [a ~ 2MeV/c^). 
The A(1520) parameters obtained from the fit are consistent with the PDG values The 
A(1520) yield, defined as the signal p.d.f. integral over the 1.48-1.56 GeV/c^ mass interval 
(2.5r), is 15519 ±412 events. 

The pKs mass spectrum is fitted to a sum of a G"*" signal p.d.f. and a third order 
polynomial. The G"*" signal shape can be rather complicated because of possible rescattering 
of particles inside nuclei In order to compare our result with other experiments we 

assume that the signal is narrow and its shape is determined by the detector resolution 
(~ 2MeV/c2). For m = 1540MeV/c2 the fit result is 29 ± 65 events. Using the Feldman- 
Cousins method of upper limit evaluation we obtain < 94 events at the 90% CL. 
We set an upper limit on the ratio of G"*" to A(1520) yields corrected for the efficiency and 
branching fractions: 

A^observed(G + ) e{pK-) S(A(1520) ^ pK') _ ^ (gQO/^ CL). 



A^observcd(A(1520)) e{pKs) B{Q+ ^ pKs) ■ B{Ks ^ tt+tt- 



It is assumed that S(G+ pKs) = 25%. We take ;B(A(1520) ^ pR-) = i-B(A(1520) ^ 
NK) = i(45 ± 1)% [9]. The ratio of efficiencies for G+ pKs and A(1520) ^ pR- of 
37% is obtained from the Monte Carlo (MC) simulation. Our limit is much smaller than 
the results reported by many experiments which observe G"*". For example it is two orders 
of magnitude smaller than the value reported by the HERMES Collaboration We do 
not know any physical explanation for such a difference. 

The momentum spectrum of the produced A(1520) is shown in Fig. El (left). This spectrum 
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FIG. 3: Left: momentum spectrum of the A(1520). Right: momentum spectrum of K (points 
with error bars) and (solid histogram). 



is obtained from fitting M{pK~) in momentum bins and correcting for the efficiency obtained 
from MC. The K~ should have a 440MeV/c momentum to produce A(1520) on a proton 
at rest. Even in the presence of Fermi motion with a typical momentum of 150MeV/c, 
A(1520) produced in the formation channel should be contained in the first momentum bin, 
0.4 to 0.6GeV/c^. Therefore most of the A(1520) are produced in the production channel. 
The projectiles that can produce A(1520) are K~ , Ks, Ki, A. The momentum spectra 
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of K~ and are given in Fig. El (right). The spectra are corrected for efficiency and 
for contamination from other particle species. It is not hkely that A(1520) production 
is dominated by interactions induced by A projectiles, because the A (1520) momentum 
spectrum is too soft. Even at the threshold of the AiV A(1520)p reaction the A(1520) 
momentum is ~ 1.1 GeV/c. 

To demonstrate that non-strange particles do not produce A(1520) we study the pK~ 
vertices accompanied by a tag. The distance from the pK~ vertex to the nearest is 
plotted in Fig. I^as a dashed histogram. For comparison the distance to any track is plotted 
as a solid histogram. The peak at zero corresponds to the vertices with additional tracks. 




Distance, cm 

FIG. 4: Distance from pK~ secondary vertex to the nearest track (solid histogram) and to the 
nearest (dashed histogram). 

The much smaller peak at zero for the tagged vertices leads us to the conclusion that 
most A(1520) are produced by strange projectiles. 



SEARCH FOR PENTAQUARKS IN B MESON DECAYS 

In this analysis we search for 0+ and 0*++ (an isovector pentaquark predicted in some 
models [l^]) in the decays — >■ Q^p followed by G"*" — *■ pKs, and Q*^~^p followed by 

0*"*"+ —y pK~^, respectively (inclusion of charge conjugated modes is implied throughout this 
section). We also search for 6° in the decay 5° Q^pn^ followed by 0° D^*^^p, and 
6*+ (the charmed analogue of 6*++) in the decay 5° 6*+p followed by 6*+ D^p. We 
reconstruct D mesons in the decay modes D*^ D^ir^, K~tt^ and D~ — >■ K't^^t:^ . 

The dominant background arises from the continuum e^e~ — > qq process. It is suppressed 
using event shape variables (the continuum events are jet-like, while the BB events are 
spherically symmetric). 

The B decays are identified by their CM energy difference, /S.E = - Ei) — -Ebeam, and 
the beam constrained mass, Mbc = \/-E'beam ~ ^^iPiYi where -Ebeam is the beam energy 
and Pi and Ei are the momenta and energies of the decay products of the B meson in the 
CM frame. The distribution (with Mbc > 5.27GeV/c2) and Mbc distribution (with 
|Ai?| < 0.05 GeV/c^) for the B^ ppKs and B^ ppK^ decays are shown in Fig. El The 
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signal yields are extracted by performing unbinned maximum likelihood fits to the sum of 
signal and background distributions in the two dimensional (Mbc, AE) space. The signal 
distributions are determined from MC, whereas the background distributions are determined 
from the AE and Mbc sideband data samples. The fits give 28.6^*^51 and 216.5]'^}g;g signal 
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FIG. 5: AE and M^c distributions for ppKg (top), and ppK^ (bottom) modes. The curves 
represent the fit results. 



yields for the ppKs and ppK^ modes, respectively. For the region 1.53 GeV/c^ < MpXg < 
1.55 GeV/c^, corresponding to the reported mass, we find no signal. Since there is only a 
theoretical conjecture for the 6*+"'', we check the 1.6 GeV/c^ < MpK+ < 1-8 GeV/c^ region 
and find no signal. Assuming both states are narrow, we set the upper limits 



B{B^ Q+p) X i3(0+ ^ pKs) 
B{B^ ppKs) 

B{B+ Q*++p) X g(e*++ ^ pK+) 
B{B+ ppK+) 



< 22% (90% CL) and 

< 2% (90% CL). 
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FIG. 6: AE and M{D'^*^p) distributions for D-ppTT+ (left), D*-ppTT+ (middle) 

and B^ — > D^pp (right) decays. The hatched histogram in AE distributions corresponds to the 
D meson sidebands, while the hatched histogram in M{D^*^p) distributions corresponds to the 
AE sidebands (shown with vertical lines on the AE plots). The vertical line in the M{D^*^p) 
distributions shows the HI mass, 3.099 GeV/c^. 



The AE and corresponding M(D^*^p) plots for the decays — > D'ppn^, — > 
D*^pp7r~^ and B^ D^pp are shown in Fig. El From the fit to AE spectra the numbers of 
reconstructed B decays are 303 ±21, 60 ± 8 and 66 ± 9 for the three modes, respectively. 
No signal of G° or G*+ is found in the M{D^*^p) spectra. We set the following upper limits 
on the fractions of the final state proceeding via G° and G*^: 
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< 1.2% (90% CL), 

< 11% (90% CL), 

< 5.9% (90% CL). 

We assume here that the charmed pentaquark mass is 3.099 GeV/c^ and that the signal p.d.f. 
is determined by the detector resolution (~ 3.5MeV/c^) j^. Our limits can be compared 
with the HI claim that about 1% of D*^ mesons originate from G|! decays. The branching 
fraction for G° decays to D mesons is expected to be even larger because of the larger phase 
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space. 
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